The Annals of Applied Probability 
2009, Vol. 19, No. 2, 556-584 
DOI: 10.1214/08-AAP552 

© Institute of Mathematical Statistics, 2009 



PORTFOLIO CHOICE WITH JUMPS: 
A CLOSED-FORM SOLUTION 

By Yacine Ait-Sahalia/ Julio Cacho-Diaz and T. R. Hurd^ 

Princeton University, Princeton University and McMaster University 

We analyze the consumption-portfolio selection problem of an in- 
vestor facing both Brownian and jump risks. We bring new tools, in 
the form of orthogonal decompositions, to bear on the problem in or- 
der to determine the optimal portfolio in closed form. We show that 
the optimal policy is for the investor to focus on controlling his ex- 
posure to the jump risk, while exploiting differences in the Brownian 
risk of the asset returns that lies in the orthogonal space. 

1. Introduction. Economists have long been aware of the potential ben- 
efits of international diversification, while at the same time noting that 
the portfolios held by actual investors typically suffer from a home bias 
effect, meaning that those portfolios tend to be less diversified interna- 
tionally than would be optimal according to portfolio choice theory [see, 
e.g., Solnik (1974) and Grauer and Hakansson (1987)]. One possible ex- 
planation is due to the risk of contagion across markets in times of cri- 
sis, which is a well-documented phenomenon [see, e.g., Longin and Solnik 
(2001), Ang and Chen (2002), Bae, Karolyi and Stulz (2003) and 
Hartmann, Straetmans and de Vries (2004)]. 

To generate this contagion, jumps of correlated sign across markets are a 
natural source of asymmetric correlation. Namely, when a downward jump 
occurs, negative returns tend to be experienced simultaneously across most 
markets, which then results in a high positive correlation in bear markets. 
When no jump occurs, the only source of correlation is that generated by 
the driving Brownian motions and will typically be much lower. 

Studying dynamic portfolio choice has a long history, going back to Mer- 
ton (1971). Subsequent papers have considered the portfolio problem, either 
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in the one-period Markowitz setting or in the more complex dynamic Mer- 
ton setting, when asset returns are generated by jump processes, including 
Poisson processes, stable processes or more general Levy processes.'^ But so 
far, when jumps are included, the determination of an optimal portfolio has 
not been amenable to a closed-form solution, and this is a long-standing 
open problem in continuous-time finance. As a result, with n assets, one 
must solve numerically an n-dimensional nonlinear equation. This is diffi- 
cult, if not impossible, to do using existing methodologies: typical numerical 
examples consider the case n = 1 where a single risky asset is present. But 
with more efficient global markets, capital flows and a considerably larger 
number of potential assets to invest in, investors have more investment op- 
portunities than ever before. In practice, we would certainly like to be able 
to solve models with hundreds if not thousands of assets. 

This paper's main contribution is to make modeling choices and bring 
powerful yet simple mathematical tools, in the form of orthogonal decompo- 
sitions, to bear on the problem: we will obtain the solution in closed form, 
thereby removing any restriction on the number of assets involved, as well 
as providing new insights and intuition into the structure of an optimal 
portfolio when jump risk is present. 



^Early papers include Aase (1984), Jeanblanc-Picque and Pontier (1990) and Shirakawa 
(1990). More recently, see Han and Rachev (2000) and Ortobelli et al. (2003) for a study 
of the Markowitz one-period mean-variance problem when asset returns follow a stable- 
Paretian distribution; Kallsen (2000) for a study of the continuous-time utility maximiza- 
tion in a market where risky security prices follow Levy processes, and a solution (up to 
integration) for power, logarithmic and exponential utility using the duality or martin- 
gale approach; Choulli and Hurd (2001) give solutions up to constants of the primal and 
dual Merton portfolio optimization problem for the exponential, power and logarithmic 
utility functions when a risk-free asset and an exponential Levy stock are the invest- 
ment assets; Liu, Longstaff and Pan (2003) study the implications of jumps in both prices 
and volatility on investment strategies when a risk-free asset and a stochastic-volatility 
jump-diffusion stock are the available investment opportunities; Emmer and Kliippelberg 
(2004) study a continuous-time mean-variance problem with multiple assets; Madan (2004) 
derives the equilibrium prices in an economy with single period returns driven by expo- 
sure to explicit non-Gaussian systematic factors plus Gaussian idiosyncratic components. 
Cvitanic, Polimenis and Zapatero (2008) propose a model where the asset returns have 
higher moments due to jumps and study the sensitivity of the investment in the risky 
asset to the higher moments, as well as the resulting utility loss from ignoring the pres- 
ence of higher moments. Das and Uppal (2004) evaluate the effect on portfolio choice of 
systemic risk, defined as the risk from infrequent events that are highly correlated across 
a large number of assets. They find that systemic risk reduces the gains from diversifying 
across a range of assets, and makes leveraged portfolios more susceptible to large losses. 
Aug and Bekaert (2002) consider a two-regime model in a discrete-time setting, one with 
low correlations and low volatilities, and one with higher correlations, higher volatilities 
and lower conditional means. 



PORTFOLIO CHOICE WITH JUMPS 



3 



Specifically, we model asset returns as following exponential Levy pro- 
cesses, which are a natural generalization of Merton's geometric Brownian 
motion model. With n assets, the space of returns is M". We show that the 
jump risk will occur in a well-defined portion of M", say V, and we then 
decompose into V and the orthogonal space, V-^. Brownian risk, by con- 
trast, occurs throughout M". By adopting a factor structure on the Brownian 
variance-covariance of returns, we obtain a matching orthogonal decomposi- 
tion of the Brownian risk which leads to a clean split of that Brownian risk 
into a portion occurring in V, where it adds up to the jump risk, and a por- 
tion in V^, where it is the sole source of risk. We are then able to distinguish 
between the optimal portfolio positions in the space V spanned by the jump 
risk (which the investor will attempt to limit) and those in the orthogonal 
space y-*" (where the investor will seek to exploit the opportunities arising 
from the traditional risk-return trade-off). 

In our model, the structure of the Brownian variance-covariance matrix 
is taken to reflect the existence of one or more economic factors or sectors, 
each sector comprising a large number of related companies or countries. 
We start with the case where there is a single economic sector and jump, 
and then consider the more general case where the economy consists of m 
sectors or regions of the world, each consisting of k firms or countries. There, 
we allow for multiple jumps, potentially as many jump terms as there are 
sectors; those jumps can affect only their sector, or all or just some of the 
sectors, and do so to different degrees. The number of sectors can be as 
large as desired, providing a fair amount of flexibility in modeling assets 
with different characteristics and sources of jump risk. 

Our closed-form solution provides new insights into the structure of the 
optimal portfolio in the presence of jumps. We show that the optimal invest- 
ment policy can be summarized by a three-fund separation theorem, in the 
case of a single jump term. In general, if there exists enough cross-sectional 
variability in the expected excess returns, then the expected return and 
volatility of the portfolio value grow linearly in the number of assets. This 
happens because of exposure to the risky assets that the investor acquires 
in the space V^. But the optimal policy in the subspace where jump risk 
lies, V, is to control the portfolio's exposure to jumps by keeping it bounded 
as the number of assets grows. As a result, the overall exposure to jumps is 
dwarfed by the exposure to diffusive risk asymptotically in the number of 
assets n. Indeed, the additional investments in the risky assets are entirely 
in the direction V-^ that is orthogonal to the jump risk: they are all achieved 
with zero net additional exposure to the jump risk. 

In other words, the optimal investment policy is to set the overall exposure 
to jump risk in V to the desired level, and then exploit, in the orthogonal 
space V'^, any perceived differences in expected returns and Brownian vari- 
ances and covariances. But if the expected excess returns have little variabil- 
ity in the orthogonal space V-^, the opportunities for diversiflcation effects 
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are weak since controlling the exposure to jumps trumps other concerns 
(including the usual Brownian diversification policy). 

The rest of the paper is organized as follows. In Section 2, we present 
our model of asset returns, and examine the investor's portfolio selection 
problem. In Section 3, we start for expository reasons with a simple case: a 
one-sector economy where the n risky assets have the same jump size and 
introduce the mathematical tools we will use to derive the optimal portfo- 
lio weights in closed form. In Section 4, we provide a comparative statics 
analysis of the dependence of the optimal portfolio weights on the arrival 
intensity of the jumps, their magnitude, and the degree of risk aversion of 
the investor. In Section 5, we extend the model to an m-sector economy 
where sectors have their own sources of jump risk and solve the optimal 
portfolio problem in that case, again in closed form. In Section 6, we extend 
the model to other utility functions. Conclusions are in Section 7. 

2. The portfolio selection model. 

2.1. Asset return dynamics. Consider Merton's problem of maximizing 
the infinite-horizon expected utility of consumption by investing in a set of 
risky assets and a riskless asset. That is, the investor selects the amounts 
to be held in the n risky assets and the riskless asset at times t £ [0,oo), 
as well as his consumption path. The available investment opportunities 
consist of a riskless asset with price So^t and n risky assets with prices 

St = ['S'l,t, . . . , Sn,t]'- 

For ease of exposition, we consider first the case of a single source of jump 
risk. In Section 5 below, we will generalize the model to include multiple 
jumps terms. Asset prices follow the exponential Levy dynamics 



with a constant rate of interest r > 0. = . . . , Wn^t]' is an n-dimensional 

standard Brownian motion. 1^ is a Levy pure jump process with Levy mea- 
sure \v[dz), where A > is a fixed parameter and the measure v satisfies 
/]gmin(l, < oo, so the jumps have finite variation. For any mea- 

surable subset A of the real line, Xi'{A) = A u{dz) measures the (possibly 
infinite) expected number of jumps, per unit of time, whose size belongs to 
A. If ^{A) = oo, it is necessarily because of small jumps, since i^iA) < oo as 
long as A does not contain 0. In other words, the expected number of jumps 
per unit of time of magnitude greater than any fixed e > is always finite. 
But there may be an infinite number of jumps of size less than e. 



(1) 



rdt 





i = 1, . . . ,n 
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The economy-wide jump amplitude Yt is scaled on an asset-by-asset basis 
by the scaling factor Jj. We assume that the support of v is such that 
Si remains positive, consistent with the limited liability provision. We also 
assume that the jump process Y and the individual Brownian motions in 
W are mutually independent. 

In the special case where y is a compound Poisson process, there are a 
finite number of jumps of all sizes and Yt = J2n=i -^"i where Nf is a scalar 
Poisson process with constant intensity parameter A > 0, and the Z^s are 
i.i.d. random jump amplitudes with law v{dz), that are independent of Nt- 
Then S follows a jump-diffusion. 

In the absence of a jump term altogether, we obtain the usual geometric 
Brownian motion dynamics for S. 

The quantities and Jj are constant parameters. We write R = 

. . . , i?n]', J = [^1, • ■ • , Jn]' ■, and S = cTfj', where 

\ ■■. \ 

We assume that S is a nonsingular matrix. 

2.2. Wealth dynamics and expected utility. Let wo,t denote the percent- 
age of wealth (or portfolio weight) invested at time t in the riskless asset 
and ujt = [uji^t, ■ ■ ■ ■,^n,t\' denote the vector of portfolio weights in each of 
the n risky assets, assumed to be adapted predictable cadlag processes. The 
portfolio weights satisfy 

n 

(4) (Jo.t + = 1- 

i=l 

The investor consumes continuously at the rate Ct at time t. In the ab- 
sence of any income derived outside his investments in these assets, the 
investor's wealth, starting with the initial endowment Xq, follows the dy- 
namics 



(5) 



dXt = -Ctdt + coo.Xt^ + 

= {rXt + uj'tKXt -Ct)dt + Xtuj[a dWt + Xtu;[J dYf 



The investor's problem at time t is then to pick the consumption and port- 
folio weight processes {Cs,uJs}t<s<oo which maximize the infinite horizon, 
discounted at rate /3, expected utility of consumption 



(6) V{Xt,t)= max Et 



OO 

e"l^'U{Cs)ds 
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subject to the dynamics of his wealth (5), and with Xt given. We wiU consider 
in detail in the rest of the paper the case where the investor has power utility, 
and briefly in Section 6 the cases of exponential and log utilities, respectively. 

Using stochastic dynamic programming and the appropriate form of Ito's 
lemma for semimartingale processes, the Hamilton- J acobi-Bellman equa- 
tion characterizing the optimal solution to the investor's problem is 

dv{Xt,t) 







max 

{Ct,wt} 



dt 



(7) 



, l d^V{Xt,t) ^, 
+ 2 dX^ 

+ xj [ViXt + Xtuj[Jz,t)-V{Xt,t)]i^{dz) 

with the transversality condition limt_»oo E[V{Xt, t)] = [see Merton (1969)]. 

The standard time-homogeneity argument for infinite-horizon problems 
gives that 



e^V(Xi,t) 



max Et 



max Et 

{Ct+u,i^t+u;0<u<oo} 



max Eo 

{Cu,Wi,;0<u<oo} 



-f3{s-t) 



U{Cs)ds 



e-f^U{Ct+u)du 
e-^'^U{Cu)du 



^L{Xt), 

which is independent of time. The third equality in this argument makes use 
of the fact that the optimal control is Markov. Thus V{Xt,t) = e~^^L{Xt) 
and (7) reduces to the following equation for the time-homogeneous value 
function L: 



0: 



max 



{u{Ct) - f3L{Xt) + ^^^{rXt + Jj\,Xt - Ct 



(8) 



+ \j [L{Xt + XtJt3z) - L{Xt)]u{dz 

with the transversality condition 
(9) lim E[e-f^*L{Xt)\ = 0. 

t— >oo 
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The maximization problem in (8) separates into one for Ct, with first- 
order condition 

dUjCt) ^ dLjXt) 
dC dX 

and one for ujt. 

\ dL{Xt) , l d^L{Xt) ^^ , 

(10) 

+ A j[L{Xt + XtJtiz)-L{Xt)]v{dz)\. 
Given wealth Xf, the optimal consumption choice is therefore 

(11) c; 

In order to determine the optimal portfolio weights, wealth and value func- 
tion, we need to be more specific about the utility function U. 

2.3. Power utility. Consider an investor with power utility, U{c) = c^~^ /(I - 
7) for c> and U{c) = -00 for c < with CRRA coefficient 7 G (0, 1) U 
(1,00). (In Section 6, we briefly treat the exponential and log utility cases.) 
We will look for a solution to (8) in the form 

(12) L(x) = K-V-7(l-7), 
where X is a constant, so that 

(13) ^ = (1 - l)L{x)/x, = -7(1 - l)L{x)/x\ 
Then (8) reduces to 

0= max (u{Ct)-(3L{Xt) + {l-j)L{Xt){Lo't'R + r) 
{Ct,wt} V 



'dU' 


-1 


'dL{Xt)\ 




{ 


V dX ) 



(14) 



^^il-j)LiXtHJ:u;t 



+ Xj[{l + u;[J z)^-"'L{Xt) - L{XtMdz 



that is. 



(15) 



0= min (--Jli^^ + -±--(r + u;iR) + ^ 

{c„u.^A (i-7)^(^t) (1-7) * Xt 



+ \Wt^u;t - J [(1 + co'^Jz)'-' - iMdz)^ , 

after division by —(1 — ^)L{Xt) < 0, so that max becomes min. 
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2.4. Optimal policies. The optimal policy for the portfolio weights ujt is 

(16) ujI = argmin5(tjj), 

where the functions 

(17) g{Lo)=-u;'R+^u;'T,uj + \'ip{u;'J) 
and 

(18) iPico'J) = -^^^ J [(1 + ^'J^)'-^ - iHdz) 
are both convex. 

Since (7,R, S,J) are constant, the objective function g is time indepen- 
dent, so it is clear that any optimal solution will be time independent. Fur- 
thermore, the objective function is state independent, so any optimal solu- 
tion will also be state independent. In other words, any optimal u>^ will be 
a constant u>* independent of time and state. Finally, the objective function 
g is strictly convex, goes to +oo in all directions, and hence always has a 
unique minimizer. In the pure diffusive case, A = and we obtain of course 
the familiar Merton solution 

(19) u;* = -S~^R. 

7 

As to the optimal consumption policy, with [dU / dC]~^ (y) = y~^^"' and 
dL{x)/dx =K~^x~^ in (11), we obtain 

(20) Q = KXt 

provided that Xt > 0. Next, we evaluate (15) at the optimal policies (C(*,u;*) 
to identify the constant K: 

(21) K = ^- iiull [^*'R + r] + -{l- 7)a;*'Sw* + ~ i^{oj*' 3). 

7 7 2 7 

The constant K will be fully determined once we have solved below for the 
optimal portfolio weights, uj* . 

Finally, we have to check that the transversality condition is satisfied. By 
plugging the optimizers X* and into (6), and then taking expectations, 
one finds 



E[V{X;,t)]=E 



oo 

e~ 



'^'U{C*)ds 



Now e-^'U{Ct)=KV{Xl,s) from which it follows that E[V{Xt,t)] solves 
df /dt = —Kf and hence decays exponentially to zero as t ^ oo, for any 
K>0. 

By plugging the constants K,uj* into (5), we see that an investor with 
power utility who selects this optimal portfolio will achieve a wealth process 
X^ which follows a geometric Levy process. 
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3. Optimal portfolio in a one-sector economy with homogenous jumps. 

3.1. The orthogonal decomposition. To begin, we consider the simplest 
possible case, where the n risky assets have identical jump size characteristics 



where J is a scalar, and 1 is the n-vector 1 = [1, . . . , 1]'. The key to finding 
the solution in closed form is to pick a basis of the space of returns, M", 
that isolates the direction of the jump vector. In this simple case, jumps are 
parallel to the vector 1, so we are led to consider the orthogonal decom- 
position = V(B V-^ where V is the span of 1 and V-^ is the orthogonal 
hyperplane. 

By construction, all jump risk is contained in V. However, Brownian risk 
occurs throughout M". We will then have to figure out which part of the 
Brownian risk takes place in V (call it S) and which part takes place in 
(call it S^). Indeed, the key to going further, and in particular to finding 
solutions in closed form, is to make a further assumption about the relation 
between the Brownian risk and the subspace V. We therefore make the 
following assumption on S: 

Assumption 1 (Orthogonal decomposition assumption). We assume that 
V and are invariant subspaces of the matrix S. In other words, we as- 
sume that S = S -|- S"*- maps V onto V and onto V-^. 

This is where assuming a factor structure for the Brownian variance- 
covariance matrix of returns S will be useful, as it will make the decom- 
position S = S -|- explicit, as we will see below. Throughout the paper, 
including the multisector case of Section 5 where the jump and factor struc- 
tures are much richer, we will always decompose the returns space M" into 
the direction of the jump vector (s), where both Brownian and jump risks 
occur (for now V is simply 1) and the orthogonal space [V-^, where only 
Brownian risk occurs). This decomposition is illustrated in Figure 1 in the 
simple case we are starting with. 

As will become clear below, this decomposition leads to a separation of 
the problem into two subproblems: one in V, and one in V-^. Because V-^ 
contains only Brownian risk, the solution in that space will be the Merton 
no-jump solution, except that it will be based not on the full S, but just on 
the fraction S"*" of Brownian risk that lies in that space. That part of the 
solution will be explicit, as the Merton no-jump solution always is. In V, we 
will have to deal with both sources of risk. 

But we will already have achieved dimension-reduction: the solution for 
the portfolio weights in the large-dimensional space, V-^, is known, and 
the remaining unknowns in V are low-dimensional (one-dimensional here in 




J = Jl 



10 



Y. AIT-SAHALIA, J. CACHO-DIAZ AND T. R. HURD 



fact). With more structure on the dynamics generating the asset returns, 
namely specific assumptions on the jump measure we will show how that 
last unknown (the optimal portfolio policy in V) can also be obtained in 
closed form. 

So, the key is to work with the orthogonal decomposition of M" that is 
suggested by the jumps. With that in mind, we therefore decompose the 
vector of excess returns as follows: 

(23) R = iil + R-^, 

where ^ is a scalar and is an ?i-vector orthogonal to 1. 
As for S, we assume for now the one- factor structure 

/ 1 

(24) E = 



P 



P 
1 



V--- P 

where > is the variance of the returns generated by the diffusive risk, 
and — l/(n — 1)</3<1 is their common correlation coefficient. 
The E matrix decomposes on = V ® as follows: 



(25) 



Ki-ll' + ft;2( I - -11' 
n \ n 



where I denotes the n x n identity matrix and 
(26) Ki=v'^ + v^{n-l)p, 



V: defined j 
by the direction 
of the jump 
vector(s> 



Jump risk: J 
Brownlan risk: E 



( Brownian risk only; .) 

V ^ : defined by orthogonality 
to the jump vi;cU)r(s) 



Fig. 1. Orthogonal decomposition of the returns space R" = V (B V where V is the 
span of the jump vector (s) and is the orthogonal hyperplane. By construction, only V 
contains jump risk. However, Brownian risk is contained throughout R" ; V contains the 
part E, and the part E^. 
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(27) K2 = 7;2(l-p) 

are the two distinct eigenvalues of S, ki with multiphcity 1 and eigenvector 
1 and K2 with multiphcity n — 1 . 

We then search for the optimal portfolio vector uj in the form 

(28) uj = u;l+uj^, 

where O is scalar and uj^ is an n-vector orthogonal to 1. 

3.2. Optimal portfolio separation. Given the decomposition (25), we see 
from (16)-(17) that the optimal Co* and uj-^* must satisfy 

{uj'^*,u!*) = arg min l—nutR-l — jnuj'^Ki + \il}{nujj) 
{w-'-.d)} t 2 



(29) 



And we now see that this separates into two distinct optimization prob- 
lems, one for iv^ and one for the scalar ud: 

: argmin{(7"'"(u;"'")}, 

(30) 




arg min{(jr(Lj)}, 



where 



(31) g^{oj^) = - i£j3i + \lK2UJ^'uJ^ 

Part of return in ^ , r d ■ ■ i ■ i/ i 

Part oi Brownian risk m V-^ 

(32) g{Lj) = — nOR + ^jnO'^Ki + Xil){nu)J). 

Part of return in V p^^^ ^^^^^^ZT^isk in V Jump risk 

In y-*-, things are simple because there is no jump risk by construction. 
The first-order condition for minimizing (31) is 

-R^ 7K2u;-^* = 0, 

whose solution is 

(33) u;^* = — R^ = —jj^ -R^. 

This is nothing else than the Merton no-jump solution, but in restriction 
to the space V'^. In other words, uj-^* is the solution in the no-jump space 
for the Merton problem with Brownian risk S"*- and the vector of excess 
returns R-*-. The investor is going after excess returns, subject to the usual 
continuous-risk provisions. 
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To solve (30) for the last remaining unknown lu* , we have to contend with 
both Brownian and jump risks. Below, we will show how to determine u)* in 
closed form under further assumptions on the distribution of the jumps and 
the investor's utility function. But for now, we see that the optimal portfolio 
choice is characterized by 

(34) u;*=Co*l + ^ -R^. 

3.3. Large n asymptotics. As far as the optimal solution for ui is con- 
cerned, with the change of variable zun = nu), we see that 

(35) tx7* = argmin|-TO„^ + l-jw^Ki/n + Xip{wnJ)\ ■ 

Letting n — > oo, we have that Ki/n ^ v'^p and so tz7* converges to a finite 
constant roj^ given by 

(36) CT^ = argmin|-tz7oo-R+ ^7^L'^^P+ -^V'(^^oo^))- 

This convergence is illustrated in Figure 3. 

An investor who selects this optimal portfolio will achieve a wealth pro- 
cess which follows a geometric Levy process with characteristic triple 
{b,c,f) such that the drift b and volatility c are 0(n), as long as there 
is sufficient cross-sectional dispersion of excess returns, in the sense that 
||R-'-|P = R-'-'R-'- = 0(n), while the Levy jump measure / remains 0(1) as 
n — > oo. 

This means that expected excess returns R"*" lead the optimal portfolio 
value to increasing expected returns b and variance c, both growing linearly 
in the number of assets. On the other hand, as n grows, the exposure to con- 
tagion jumps remains bounded, and is dwarfed by the exposure to diffusive 
risk. Indeed, the investment in the risky assets due to the expected return 
compensation R^ is entirely in the direction of u>^, which is orthogonal to J. 
Since all the assets have the same response to the jump risk under the simple 
model (22), this is all done in a long-short manner: the sum of the wealth 
percentages invested in is 0. So the additional amounts invested in the 
risky assets are all achieved without increasing the net exposure to the jump 
risk. In fact, since ro* — > tu^ finite, we have asymptotically (D* ~ vj*^/n so 
the investor exploits an increase in the number of assets to reduce, ceteris 
paribus, his net exposure to the risky assets. 



3.4. A three-fund separation theorem. It is clear from the analysis above 
that there exist three funds such that an investor should be indifferent be- 
tween the original n assets (plus the riskless asset) or just the three funds; 
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the proportions of each fund invested in the individual assets depend only 
on the asset return characteristics and not on investor preferences; and the 
investor's optimal demands for the three funds depend on his preferences. 
By contrast, both in the static one-period Markowitz model and in the dy- 
namic Merton model, one needs two funds (the market portfolio and the 
riskless asset) to replicate the optimal asset allocation. 

In our simple one-sector model with homogenous jump sizes, the three 
funds are: 

• A long-short fund LS which holds proportions Si = R-*- of the n 
risky assets. This fund is long-short because it contains no net exposure 
to the risky assets, since the sum of the weights invested in the various 
assets by fund LS is zero. Said differently, any investment in fund LS 
generates no new exposure to jump risk. 

• An equally weighted fund EQ holding a portfolio ^2 = of the n risky 
assets. Investment in this fund generates net exposure to jump risk. How 
much of that risk the investor is willing to take is controlled by the quan- 
tity Co*. 

• A risk-free fund RF holding the riskless asset only. 

Funds LS and RF provide an efficient diffusion risk-return optimization 
(for the part of the Brownian risk orthogonal to 1, that is, the part in V-^), 
while fund EQ allows investors to hedge against the jump risk, plus the part 
of the Brownian risk in the direction of 1, that is, the part in V. Replicating 
the investor's optimal positions cannot be achieved by lumping together LS 
and EQ into a single market portfolio. 

The fact that the "jump-hedging" fund is perfectly equally weighted (or 
equivalently the fact that the orthogonal fund is perfectly long-short) is of 
course due to the fact that in the model so far all assets are symmetrically 
exposed to the jump risk. The fact that we need three funds is due to the 
fact that so far we have only one source of jump risk. We will generalize all 
of this below in Section 5. 

3.5. Fully explicit portfolio weights. Some special cases lead to a closed- 
form solution for the last remaining constant, tu*, yielding fully closed- form 
solutions for the optimal portfolio weights. For this, we need to specify the 
utility function and the Levy measure i'{dz) driving the common jumps; 
then we can compute the integral in (36). We provide a few examples. 

Consider the case of a power utility investor with CRRA coefficient 7 = 2, 
the Levy measure satisfying a power law. 



(37) 




if zG(0,l], 
ifzG[-l,0), 
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and J G (—1,1). We define A4- > as the intensity of positive jumps and 
A_ > as tfie intensity of negative jumps. Equation (36) specializes to 

(38) w*^ = arg min /oo {w) , 

■OJ 

where 

(39) /oo(ro) = -wR + w^v^p + X+ I [(1 + wJzy^ - 1] dz/z 

Jo 

(40) ~^~J iC^ + ^J^y^ -'^]dz/z 

= —wR + ■nP'v^ p — A+ log(l + wJ) — A_ log(l — wJ\ 
The first-order condition (FOC) for vj is given by 

(41) - + Ivjv'^p - A+ J(l + vjjy^ + A_ J(l - wjy^ = 0. 
The optimal solution must satisfy 

(42) \JwU<^\ 

otherwise, there is a positive probability of wealth Xt becoming negative, 
which is inadmissible in the power utility case. The asymptotic solution to 
(41) is the unique root w*^ satisfying the solvency constraint (42), and that 
solution is given by 

* R 



6v^p 



/4/ i?2 A+ + A_ 1 
^ V 3 \l2{v^pf ^ 2v^p ^ J2 

1 // / R^ A+ + A- 



(43) 



xcosl-arccos^^^^^^^+ ^^^^ 

_ 3(A+ - A,) _ J_ 
2JR J2 



'1/ i?2 A+ + A_ 1 



3\12{v^p)^ 2v^p J2 



4 
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3.5.1. Negative jumps. Let us now consider the case J < 0, A_ = and 
A = A+ > in order to capture the downward risk inherent in the types of 
jumps we are concerned about. Figure 2 plots function of J and 

A. Visual inspection of (43) reveals that J and A do not have a symmetric 
effect on the optimal portfolio weights. In the exact small-sample case, the 
optimal solution (43) becomes 



Figure 3 plots the objective function, /^(ro) = g{w/n), and shows its con- 
vergence to /oo('C^) as n^oo, along with argmin/„(tA7) = zu*, converging 
to ro^. It is also worth noting that 



so that the optimal investment in the risky assets is always less than what 
it would be in the absence of jumps, nR/ {2ki). This is natural since J < 0. 

3.5.2. Other examples. Other cases that lead to a closed-form solution 
for tu* include power utility with 7 = 3 and either power law jumps u[ dz) = 
\dz/z, or uniform jumps u[dz) = Xdz. In either case, the FOG is then a 
cubic equation, as (41), solvable in closed form using standard methods. 
Another case is one where the investor has log utility with jumps of a fixed 
size, v{dz) = 5{z = z)dz, for some z G [—1,1]. Then the FOC leads to a 
quadratic equation. 

4. Comparative statics. Using the explicit solutions above, we can now 
investigate how the optimal portfolio responds to different jump intensities, 
jump sizes and degrees of risk aversion. This analysis is decidedly partial- 
equilibrium, in that we do not allow for changes in the assets' expected 



(44) 




Aki 



(45) 









1 - 



1 



Fig. 2. Optimal portfolio weight uj^ as a function of J and X. 
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Fig. 3. Scalar objective function used to determine the optimal portfolio weight zu^ and 
its large asset asymptotic limit, vo^. 

returns when other parameters of the model change. To save space, we con- 
sider only the negative jump example of Section 3.5.1, where J < 0, A_ = 
and n — > oo. 

4.1. Response to jumps of different arrival intensity. We have 
(46) w*^^—Qo asA— >oo, 



The first limit means that the investor will go short to an unbounded extent 
on all the risky assets if the arrival rate of the jumps goes to infinity. This is 
to be expected, since J < and we impose no short sale constraints. Further, 
ro^ tends to — c« when A ^ oo at the following rate: 



If, on the other hand, the jumps become less and less frequent, A — > 0, 
then tz7^ tends to a finite limit driven by the diffusive characteristics of the 
assets, R/{2v'^p), which is the limit in the no-jump case, unless the jump 
size J is so large that the solvency constraint binds. This gives rise to a kink 
in the demand function. 

If the solvency constraint is not binding, then the higher the variance 
of the assets and/or the more heavily correlated they are, the smaller the 
investment in each one of them. And the higher the expected excess return 



(47) 




as A 



0. 
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of the assets the higher the amount invested. For a small perceived jump 
risk (A small), the optimal solution behaves like 

4.1.1. Jumps versus expected return trade-off. The weights are mono- 
tonic in A, with 

^J{2v'^p - JRf + 8 J(i? + JA)w2p 

If > 0, there exists a critical value A such that 

(48) > for A < A, 

(49) < for A > A. 

That is, as long as jumps do not occur too frequently (A < A), the investor 
will go long on the assets in order to capture their expected return, even 
though that involves taking on the (negative) risk of the jumps. When the 
jumps occur frequently enough (A > A), then the investor decides to forgo 
the expected return of the assets and focuses on canceling his exposure to 
the jump risk by going short these assets. 

The critical value A takes a particularly simple form. It is given by 

(50) A ^ 



J| /o zv{dz) 

Clearly, the higher R relative to |J|, the higher A. And the smaller the 
expected value of Z, the bigger A. With v{dz) = dz/ z, we get 

(51) A = ^. 

This expression can be interpreted as an analogue to the Sharpe ratio, but for 
jump risk: excess return in the direction of jumps, -R, divided by a measure 
of the jump magnitude, | J|. 

Now, if ^ < 0, then zz7^ < for every A > 0. In that case, there is no point 
in ever going long those assets since both the expected return and the jump 
components negatively impact the investor's rate of return. 

4.1.2. Flight to quality. The solution above can capture a well-documented 
empirical phenomenon. Imagine an investor currently in a normal, low jump 
risk environment, who receives information suggesting that the jumps be- 
come more likely. Starting from a situation where A < A, if the perception 
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of the jump risk increases (A | A), then the optimal poHcy for the investor is 
to flee-to-quahty, by reducing his exposure to the risky assets (ro^ | 0) and 
investing the proceeds in the riskless asset. If the perception of the jump 
risk exceeds the critical value A given in (50), then the investor should go 
even further and start short-selling the risky assets. Because the jump risk 
affects all the assets, the perception of an increase in the intensity of the 
jumps leads the investor to dump all the risky assets indiscriminately. 

4.2. Response to jumps of different magnitudes. If we now concentrate 
on the effect of an increase in the jump size instead of the jump intensity, 
then 

d^lo_ 1 ^ 2pv'' + JR 



^•^ 2 J2 V J{2v^p - JRY + SJ{R + J\)v'^p 



for J < 0. The monotonicity implies that as the jump size gets closer to zero 
(J t 0), the investor increases his holdings in the risky assets and conversely 
as Ji (-1). 

As to the sign of vj%^, we have 

(52) > for - R/\ < J < 0, 

(53) ■^lo<^ for - 1< J < -R/\, 

as long as.R/A<l. If-R/A>l, the expected return is high enough relative 
to the jump intensity that the investor will always maintain a positive in- 
vestment w*^ > in the different assets, no matter how large the jump size 
(within the constraint J > — 1 , of course) . 

4.3. Sensitivity to risk aversion. Here we consider the effect of the CRRA 
coefficient 7 on ro^. For a CRRA investor, the first-order condition of (36) 
is given by 

(54) - R + 'jwoov'^ p - X [ Jz{l + WooJz)~^i'{dz) = 0; 

Jo 

then, making use of the implicit function theorem, we get 

(55) _ rul,v^p + XJ^Jz{l + wl,Jz)-nn{l + w*^JzMdz) 

37 ^{v^p + Xj^,Pz^il + w*^Jz)-^-^i^{dz)) 

The denominator is always positive but the numerator could be negative, 
zero or positive depending on the sign of ro^. That is, 

> 0, if TO^ < 0, 



d 



(56) r^ = <'=0, ift:^7^=0, 

[ < 0, if w*^ > 0. 

This, in turn, implies that the higher the CRRA coefficient of an investor is, 
the smaller will be his ro^ in absolute value. In the limit where 7 increases 
to 00, |ct^| decreases to zero. 
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5. Optimal portfolio in a multisector economy with multiple jump risks. 

We now generalize our previous results by studying the more realistic port- 
folio selection problem in an economy composed of m sectors (or regions of 
the world), each containing k firms (or countries). We also generalize the 
model to allow for multiple sources of jump risk. The total number of assets 
available to the investor is n = mk. 

To understand the qualitative implications of the solution for portfolio 
allocation in the presence of jumps, one would presumably be primarily 
interested in the situation where m is fixed and k goes to infinity with n. 
But we provide the full solution, including the special case where all the 
assets are fundamentally different, that is, k = 1 and m = n. This provides 
a fair amount of generality to capture different situations, with different 
numbers of sectors. 

5.1. Sector jumps. As in the simpler one-factor model of Section 3, we 
will decompose the space of returns as M" = V ® V-^ where V is the span of 
the jump vector(s) and V-^ is the orthogonal hyperplane. With m sectors, it 
is now natural to assume that there is potentially a separate jump term per 
sector, which can affect not only its own sector but also the other sectors, 
and do so to different degrees. 

We are therefore replacing the dynamics of the risky asset returns in (2) 
with 

dS' ^ 

(57) -^ = {r + R,)dt + Y,ai.jdWj,t + Y.'Ji,i(^Yi,t^ i = l,...,n, 

j=i 1=1 

where each Yi^t is a Levy pure jump process with Levy measure Xiiyi{dz), 
I = 1, . . . ,m. We assume that the Levy pure jump processes and the Brow- 
nian motions are mutually independent and that the sum of the jumps has 
support on (— l,c>o) to guarantee the positivity (or limited liability) of Si. 
The constant Ji^i is asset i's scaling of sector Ts jump. Let J; = [Ji^i, . . . , Jn,i]' ■ 
Let us define the vector 1/ as the n- vector with I's placed in the rows 
corresponding to the /-block and zeros everywhere else: 

(58) li = [0,...,0,l^_^,0,...,0]', 

sector / 

where the first 1 is located in the k{l — 1) + 1 coordinate, and the last one 
in the kl + 1 coordinate. 

We assume that each jump vector is of the form 

m 

(59) Jl = ^js^ils = [ji^i, . . . ,ji^l,j2,l, ■ ■ ■ ,32,h ■ ■ ■ ,jm,l, ■ ■ ■ ,jm,l]'^ 
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meaning that firms within a given sector have the same response to the 
arrival of a jump in a given sector, that is, to a jump in Yi, but the propor- 
tional response jg^i of firms of different sectors to the arrival of a jump can 
be different, and also the proportional response of firms in a given sector to 
the arrival of jumps of different sectors can be different. Since some of the 
js,i coefficients can be zero, jumps in one-sector can affect only this sector 
or some or all of the sectors. 

So the jump vectors are unrestricted linear combinations of the m vec- 
tors {li}i 

=i,...,m- Corresponding to the above structure, we now define the 
orthogonal decomposition = V (B V-^ where V is the m-dimensional span 
of the vectors {1;};=!,...,™ that contains the jump vectors, and is the 
(n — m) -dimensional orthogonal hyperplane. By construction, contains 
no jump risk, while V contains all the jump risk. As in Section 3.1, to obtain 
meaningful results, we need to make the orthogonal decomposition assump- 
tion to be able to write S = S -|- S^. With this assumption, the picture is 
similar to Figure 1 except that V is now of dimension m. 



5.2. Multifactor Brownian risk. To capture the notion of sectors, we 
specify a block-structure for the variance-covariance matrix of returns. 



(60) 



nxn 



^2,1 



Si, 2 
Sm m— 1 



^2,m 



with within-sector blocks 



1 Pl,l 



(61) 



kxk 



■ • Pl,l 
Pl,l 1 



and across-sector blocks 



/ Pl,s PL: 



(62) 



^l,s 
kxk 



ViVs 



PL: 



V 



pi,: 



Pi,: 
Pi,: 



where 1 > pi^i > pi^s and pi^i > — l/{k — 1) and pi^s > —^/{n — 1). In an asset 
pricing framework, this corresponds to a multifactor model for the returns 
process with m common Brownian factors and n idiosyncratic Brownian 
shocks. 

Some structure on S is needed, because without an explicit decomposition 
of S there is no explicit optimal portfolio policy except to say that it is in the 
form (S-'-)~^R"'" but we would not know what S"*" is. The block-structure 
of E means that we have "within-sector homogeneity" of the assets, instead 
of the complete homogeneity in the model of Section 3. Since the number 
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of sectors can be arbitrarily large, this is not a serious restriction and one 
that would have to be assumed anyway when one estimates the model lest 
the number of parameters becomes unwieldy. 

With this structure on S, it turns out that we can find exactly which part 
of the T, matrix lines up with the jumps (in V defined by the specification 
of the m jump terms), and which part is orthogonal to them (in V'^). That 
decomposition of S is given by 

m 

(63) s=^^i^i'^+S\ 

l,s=l 



where 

(64) K«„ 



vfil + {k-l)p^), iil = s, 
kvivspi^s, if / / s. 



1;!^ is an n X n matrix with a k x k matrix of I's placed in the (/,s)-block 
and zeros everywhere else. 

The part of the S matrix orthogonal to the jumps is given by the block- 
diagonal matrix 



nxn 







(65) S 
with blocks 

[{k-l)/k -l/k 

(66) i:i-^i = vf{i-p^) 



kxk 



-l/k ■•• -l/k 
V ••• -1/^ ik-l)/k^ 



5.3. Optimal portfolio separation in a multisector economy. We decom- 
pose the vector of expected excess returns on the same basis as above, 

m 

(67) R = ^nli + R^=R + R^, 

1=1 

where R-*" is orthogonal to each 1; and has the form 

R-^ = i^i', ■ ■ • ,RmT- 

Each of the fc- vectors R^ is orthogonal to the fc-vector of I's. 

We will be looking for a vector of optimal portfolio weights in the form 



(68) to = y^CJ;!; -|- CJ"*" =UJ + U! 
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By the orthogonal decomposition assumption, S"*" is orthogonal to the space 
generated by {l/}/=i,...,m- Thus, the CRRA investor's minimization problem 
again separates as 



(69) 



where 




= arg min{ g^{uj'^)}, 
arg min{^(a))}, 



(70) (w^) = -c^^'R^ + i7a;^'S^a;^ 



and 



(71) 



with 



g{Cd) = Cd'R + ^w'Ew + ^iM'^'^i 



m 



1 = 1 /,S = 1 l = \ \ S = \ J 



m 



(72) Ux) = -^^^ J [(1 + ^^)'"" - ^Mdz). 

The key consequence of this separation is that in the space V'^ where only 
Brownian risk occurs, we again get back the Merton no-jump solution (with 
S"*" and R"*" instead of the full S and R). Indeed, the first-order condition 
for minimizing g-^{uj-^) is 

(73) E^Lj^* = -R^, 

7 



all /, leads to the optimal solution uj^* = [wf *', . . . ,0;^*']' with 



which, by the block diagonal form of S and by construction l/o; = for 

±*/ , ,±* 

1 > • • • > '^m 

1 

for I = 1, . . . ,m. Note that in general, the null space of S"*- is V and so (73) 
has a unique solution in V'^. 

In the space V, both jump and Brownian risks coexist. Minimizing g{uj) 
leads to the analogue of what happens with one-sector, as in Section 3, but 
in dimension m. Similarly its solution has a limit as k goes to infinity with n 
(the number of sectors m being fixed). With the change of variable zun = ku) 
we see that 




Return contribution (part in V) 
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(75) + o -T^^riZ^ns 

l,s=l 



Brownian risk contribution (part in V) 



Jump risk contribution (all in V) 

which, compared to (36), is an m-dimensional minimization problem, instead 
of a one-dimensional one. The convexity of the objective function implies the 
existence of a unique minimizer. Letting /c — > oo, we see that ro* converges 
to a finite limit, va^, with the limit given by the equation above with Ki^g/k 
replaced by viVsPi^s- As in the one-factor case, we can determine below ro* 
in closed form under some specific jump distributions. 

The wealth process of the optimizing investor will have geometric 
Levy dynamics. Under the natural condition that R^'R^ = 0{k), as k, the 
number of stocks per sector, increases, the optimal portfolio can achieve 
expected gains at the expense of variance which both grow approximately 
linearly with k, while keeping the exposure to the sector jumps bounded. 
This result is achieved by the investor apportioning an increasing fraction 
of assets in the subspace orthogonal to the vectors 1/. And it follows 
that a generalization of the separation theorem of Section 3.4 will also hold 
in this situation. 

5.4. Three-fund separation. More specifically, consider for simplicity a 
model in which the common jumps are identically distributed, that is, i'i{dz) = 
v[dz) and J; = J for all I = 1, . . . ,m. The objective is to minimize 

(76) g{Lo)=-uj''R+^Lj'T.uJ + Xmij{uj'J) 
with 

(77) i^iiu'J) = -JY^ I [(1 + ^'J^)'-' - mdz). 

The first-order conditions for to are given by 

(78) -R + 7Su; + mAJVj(w'J) =0. 

Define the scalar y = lo'J. After multiplying (78) by ^J'S"-*^ we have that y 
must satisfy 



(79) - ij'S-iR + y+— J'S-ij'0(y) = O; 

7 7 
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then use (78) to solve for uj: 

(80) a; = is-iR-— S"ij'0(y) 

7 7 



or 



(81) 



7 V J'S-iJ 

It is clear from the analysis above that three-fund separation continues 
to hold in the multisector case. The three funds are: 

• a mean- variance fund which holds proportions 5i = S~^(R— j,^ J) of 
the n risky assets, 

• a hedging fund holding a portfolio 82 = j,^^_ij S~"^J of the n risky assets; 
how much of this fund the investor is willing to take is controlled by the 
quantity y, 

• a risk-free fund holding the riskless asset only. 

When the common jumps of the m-region model are not identical, the 
three-fund separation property generalizes to an (m -|- 2)-fund separation 
property. 

5.5. Examples of fully closed-form portfolio weights. For the m-region 
model in which the common jumps are identically distributed, a CRRA 
investor with CRRA coefficient 7 will optimize the following objective func- 
tion: 



m m 



gni-cu/k) = -^Wniri + 

1=1 ^ l=\s=\ ^ 



(82) 



where 

(83) ^(y) = -JY^ I [(1 + y^)'~' - ^Hdz). 
The optimal solution must satisfy the solvency constraints 

m. m 

(84) Zsup rjJnsjs > - 1 and Zi^i ^ Wnsjs > -1, 

s=l s=l 
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where Zgup and Zi^{ are the upper and lower bounds of the support of v{dz). 
Otherwise, there is a positive probabihty of weahh Xt becoming negative, 
which is inadmissible in the power utility case. 
The first-order conditions for ro„ are given by 

m / m \ 

(85) - n + 7 X! + ji^^'^ X! '^nsjs = 

s=l ^ \s=l J 

for I = 1, . . . ,m. These first-order conditions form a system of m equations 
which admit a unique solution -ajn satisfying the solvency constraints. Let 

j = [jl,- ■ • Jm]', r = [ri,. . .,rm]', VOn = [Wnl, ■ . . ,Wnm]' , 

(^1,1 ••• Kl,r 
; ■•. ; 

and y = X^^Li '^nsjs = j'^n- Then (85) becomes 

(87) - r + |Kro„ + mAjV'(?/) = 0; 
multiplying (87) by ^j'K~^ we get that y must satisfy 

(88) - -j'K-V + y + !^j'K-ij^(y) = 0. 

7 7 

Once we know ij: we can solve (88) for y and then use (87) to solve for Wn'- 

(89) .^„ = V^r-^K-i^y) 

7 7 

or 

(90) ^„.*K-.r + K-'j(ziZ^!^^ 

7 V j'K 

We again allow the Levy measure to generate asymmetric positive and 
negative jumps. So, we consider a Levy measure such that 

with A+ > and A- > 0. For a reasonable value of 7, say 7 = 2, then 
mXip{y) = — A+ log(l + y) — A_ log(l — y) and (88) becomes a cubic equation 
in y: 

(92) fiy) = -A + y + B[-X+{1 + y)"! + A_(l - y)"^] = 
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Fig. 4. Bivariate objective function in a two-sector economy. 



with solvency constraint \y\ < 1, where A = -Ij'K^^r and B = ■|j'K~^j. The 
optimal solution to (88) under the solvency constraint is solvable fully in 
closed form: 




-1 ffAf ^ A^ 3B{X -A+) 
X cos - arccos — -IH h 



3 VV3V 9 2A 

(93) ^B(X 




1 / 42 



4 

+ -TT 

3 



Multiple sectors provide the investors with additional options. For instance, 
if r = — (A+ — A--)j, then the investor has no exposure to jump risk. The 
reason is that jumps in one-sector are used to offset jump risk in another 
sector. Figure 4 plots the objective function, fnij^) = gi'^/k) that we obtain 
in a two-sector economy, that is the function (82) with m = 2. 

6. Other utility functions. Before concluding, we briefly outline how 
comparable results can be derived for the other two main families of prefer- 
ences, namely exponential and log utilities. 

6.1. Exponential utility. As seen above, for an investor with CRRA util- 
ity function the optimal wealth process achieved by picking the constant 
portfolio fractions lo* is itself a one-dimensional geometric Levy process 
whose characteristic triple is {Xf_{LO*''R + r — K), Xf^uj*'aa'u!* , f) where 
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f{dy) = u{dz) with y = X^^u'Jz. In this case, the investor keeps constant 
fractions of wealth in each risky asset, and the constant remaining fraction 
1 — J2j the riskless asset. 

Now, consider an investor with exponential utility, U{C) = —- ex.p{—qC) 
with CARA coefficient q> 0. We can look for a solution to (8) in the form 

(94) L(x) = -I^e-""^, 

q 

so that 

(95) — ^ = -rgL(x), \ =r^q^L{x). 

ox ax^ 

Then (8) reduces to 

= max [u{Ct) - pL{Xt) - rqL{Xt){Xtr + XtUj[Yi - Ct) 

{Ct,u)t] L 

(96) + ]-r''q^L{Xt)X^^u:[i:uJt 



2 



+ \j [e-^^^^-^J^L(XO - L{Xt)]v{dz)Y 



that is, 



(97) 



= min (i^^ - ^ - {Xtr + X^^^R - Ct) 
{Ct,Lot}VrqL{Xt) qr ^ 

+ \rqX^u:[i:u:t + ^ / [e'^^^-^-^^ - l]v{dz)^ 

after division by qL{Xt) [note that max becomes min as a result of qL{Xt) < 
0]. 

The optimal policy of u; = XfUt is given by the objective function 

(98) minf-a;'R + -rquj'EuJ + — f [g-'^'^'-'^ - l]iy{dz 
{u;}\ 2 rqj 

and the optimal consumption choice is therefore 

(99) C; =rXt-^log{rK). 
Finally, we evaluate (97) at C* and lo* to identify K, 

(100) = - exp f 1 - ^ - qLO*'-R+ -rq'^io*'i:uj* + - / [e"'?^*'-^^ - l]u(dz 

r \ r 2 r J 

For an investor with CARA utility function the optimal wealth process X^ 
achieved by picking constant portfolio amounts is now a one dimensional 
arithmetic Levy process. The investor keeps constant amounts of wealth in 
each risky asset, and the remaining amount in the riskless asset. 
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6.2. Log utility. Finally, consider an investor with log utility, U{x) 
log(a;). We can look for a solution to (8) in the form 

(101) L{x) = K{hog{x) + K2, 

where Ki and K2 are constant, so that 

Then (8) reduces to 

= max (log(Q) - (3K^^ log(Xt) - (3K2 

{Ct,OJt} I 

(103) + K^'X^\rXt + LoiRXt - Q) 

- ^K^^uj't^uJt + AiTf 1 J log(l + uj'tJz)i^idz)y 



that is, 



min (- log(a) + (3K{^ log(Xt) + 
Ct,uit} I 



{Ct 

(104) - i^f V - A'f^w^R + K^^X^^Ct 

+ A'f^icj^Swj - XK{^ J logil+io'tJ z)i^{dz)y 
The optimal policy of ujt is given by the objective function, 

(105) min f-WjR + l-Lo[T.LOt - A / log(l + uJt3z)i^{dz)] 
{uJt} \ 2 J J 

and the optimal consumption choice is therefore 

(106) C: = KiXt. 

To identify Ki and K2, we evaluate (104) at C* and u!*, 

K2 = ^{log(/5) + ^ + 1 - + 1/ log(l +u;*'Jz)z.(dz)|, 

7. Conclusions. We have proposed a new approach to characterize in 
closed form the portfolio selection problem for an investor concerned with 
the possibility of jumps in asset returns, and who seeks to control this risk 
by diversification or other means. We extended the standard multiasset ge- 
ometric Brownian motion models to exponential Levy models through the 
inclusion of correlation effects due to jumps. 
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Our key decomposition is one of the space of returns into a space contain- 
ing both jump and Brownian risks, and one containing only Brownian risk. 
The model is quite flexible and can capture many different situations, from 
jumps of different signs to jumps of a single sign, multiple jumps that can 
affect one or some of the sectors; it is as flexible as it can get while achieving 
a closed-form solution, which is critical if one is to solve a portfolio choice 
problem in practice with a large number of assets. 

With explicit optimal portfolios policies, one can address important prac- 
tical questions. How exactly does increasing the number of available assets 
improve the investor's exposure to both diffusive and jump risk? How does 
the portfolio of an investor who fears jumps differ from the portfolio of one 
who does not? Is there a simple form for the optimal portfolio which is 
achieved asymptotically as the number of assets grows to infinity? And so 
on. 

Finally, we discuss further possible extensions and generalizations of this 
work: 

1. One could consider modeling the spectral decomposition of the S matrix 
directly, instead of parametrizing the matrix itself and then determining 
its decomposition across the two spaces V and V-^. As an empirical strat- 
egy, one could determine the number of sectors through factor analysis 
or similar techniques, in order to determine in a data-driven manner the 
shape of the S matrix. 

2. Stochastic volatility of the type considered in Liu, Longstaff and Pan 
(2003) requires solving our nonlinear equations for weight vectors stepwise 
in time, in parallel with ordinary differential equations (which themselves 
depend on the current portfolio weights). This does not appear doable in 
closed form. 

3. Portfolio restrictions such as short-selling constraints are relevant in prac- 
tice but, when generic constraints are imposed on the optimal portfolio, 
we cannot expect the dimensional reduction to be preserved or our con- 
clusions to hold. However, a utility function such as power utility which 
becomes — oo for wealth below a finite threshold, sometimes automati- 
cally implies certain constraints; it appears that in this case, much of our 
analysis remains intact. 

4. Finally, we would like to be able to capture a more subtle form of conta- 
gion, in the form not just of simultaneous jumps within or across sectors, 
as we are currently able to model, but rather in the form of a jump 
in one-sector causing an increase in the likelihood that a different jump 
will occur in another sector. To capture this effect, self-exciting jump 
processes seem a promising approach which we intend to investigate in 
future work. 



30 



Y. AIT-SAHALIA, J. CACHO-DIAZ AND T. R. HURD 



Acknowledgments. We are very grateful to Jun Liu, Francis Longstaff 
and Raman Uppal for comments. We are also grateful to seminar partici- 
pants at LBS, LSE, Stanford, Tilburg, UIC, UIUC, Wharton, the European 
Science Foundation Conference on Advanced Mathematical Methods for Fi- 
nance and the Econometric Society Australasian Meeting (E. J. Hannan 
Lecture) . 

REFERENCES 

Aase, K. K. (1984). Optimum portfolio diversification in a general continuous-time model. 

Stochastic Process. Appl. 18 81-98. MR757349 
Ang, a. and Bekaert, G. (2002). International asset allocation with regime shifts. Review 

of Financial Studies 15 1137-1187. 
Ang, a. and Chen, J. (2002). Asymmetric correlations of equity portfolios. Journal of 

Financial Economics 63 443-494. 
Bae, K.-H., Karolyi, G. a. and Stulz, R. M. (2003). A new approach to measuring 

financial contagion. Review of Financial Studies 16 717-763. 
Choulli, T. and Hurd, T. R. (2001). The role of Hellinger processes in mathematical 

finance. Entropy 3 150-161. 
CviTANiC, J., POLIMENIS, V. and Zapatero, F. (2008). Optimal portfolio allocation with 

higher moments. Annals of Finance 4 1-28. 
Das, S. and Uppal, R. (2004). Systemic risk and international portfolio choice. Journal 

of Finance 59 2809-2834. 
Emmer, S. and Kliippelberg, C. (2004). Optimal portfolios when stock prices follow an 

exponential Levy process. Finance Stock. 8 17-44. MR2022977 
Grauer, R. and Hakansson, N. (1987). Gains from international diversification: 1968- 

1985 returns on portfolios of stocks and bonds. Journal of Finance 42 721-739. 
Han, S. and Rachev, S. (2000). Portfolio management with stable distributions. Math. 

Methods Oper. Res. 51 341-352. MR1761863 
Hartmann, p., Straetmans, S. and de Vries, C. (2004). Asset market linkages in crisis 

periods. Review of Economics and Statistics 86 313-326. 
Jeanblanc-Picque, M. and Pontier, M. (1990). Optimal portfolio for a small investor in 

a market model with discontinuous prices. Appl. Math. Optim. 22 287-310. MR1068184 
Kallsen, J. (2000). Optimal portfolios for exponential Levy processes. Math. Methods 

Oper. Res. 51 357-374. MR1778648 
Liu, J., Longstaff, F. and Pan, J. (2003). Dynamic asset allocation with event risk. J. 

Finance 58 231-259. 

Longin, F. and Solnik, B. (2001). Extreme correlation of international equity markets. 
J. Finance 56 649-676. 

Madan, D. (2004). Equilibrium asset pricing with non-Gaussian factors and exponential 
utilities. Technical report, Univ. Maryland. 

Merton, R. C. (1969). Lifetime portfolio selection under uncertainty: The continuous- 
time case. Review of Economics and Statistics 51 247-257. 

Merton, R. C. (1971). Optimum consumption and portfolio rules in a continuous-time 
model. J. Econom. Theory 3 373-413. MR0456373 

Ortobelli, S., Huber, L, Rachev, S. T. and Schwartz, E. S. (2003). Portfolio choice 
theory with non-Gaussian distributed returns. In Handbook of Heavy Tailed Distribu- 
tions in Finance (S. T. Rachev, ed.) 547-594. Elsevier, Amsterdam. 

Shirakawa, H. (1990). Optimal dividend and portfolio decisions with Poisson and 
diffusion-type return process. Technical report, Tokyo Institute of Technology. 



PORTFOLIO CHOICE WITH JUMPS 



31 



SOLNIK, B. (1974). Why not diversify internationally rather than domestically? Financial 
Analysts Journal 30 48-53. 

Y. Ait-Sahalia J. Cacho-Diaz 

Department of Economics Department of Economics 

Princeton University Princeton University 

AND Princeton, New Jersey 08544-1021 

NBER USA 

Princeton, New .Jersey 08544-1021 E-mail: jcacho@princeton.edu 

USA 

E-MAIL: yacine@princeton.edu 

T. R. Kurd 

Department of Mathematics 

AND Statistics 
McMaster University 
Hamilton, Ontario L8S 4K1 
Canada 

E-MAIL: hurdt@mcmaster.ca 



